Numerical study of Kondo impurity models with strong potential scattering: 
- reverse Kondo effect and antiresonance - 
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Accurate numerical results are derived for transport properties of Kondo impurity systems with potential 
scattering and orbital degeneracy. Using the continuous-time quantum Monte Carlo (CT-QMC) method, static 
and dynamic physical quantities are derived in a wide temperature range across the Kondo temperature 7k. 
With strong potential scattering, the resistivity tends to decrease with decreasing temperature, in contrast to the 
ordinary Kondo effect. Correspondingly, the quasi-particle density of states obtains the antiresonance around 
the Fermi level. Thermopower also shows characteristic deviation from the standard Kondo behavior, while 
magnetic susceptibility follows the universal temperature dependence even with strong potential scattering. It is 
found that the t-matrix in the presence of potential scattering is not a relevant quantity for the Friedel sum rule, 
for which a proper limit of the /-electron Green's function is introduced. The optical theorem is also discussed 
in the context of Kondo impurity models with potential scattering. It is shown that optical theorem holds not 
only in the Fermi-liquid range but also for large energies, and therefore is less restrictive than the Friedel sum 
rule. 
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PACS numbers: 

I. INTRODUCTION 

A single magnetic impurity embedded into the sea of con- 
duction electrons shows Kondo effect. Although this prob- 
lem has been almost continuously studied during the last 40 
years, dynamic and magnetic properties of the Kondo and re- 
lated models still attract great interest in condensed matter 
physics. Much less attention is paid to anomalous decrease 
of resistivity with lowering temperature since such decrease 
can be caused by many different mechanisms. The first at- 
tention to this problem is traced back to the 60's, when resis- 
tivity of dilute Fe and Cu alloys in Rh_matrix revealed a new 
type of anomaly at low temperatures!! Namely, it was ob- 
served that resistivity decreases with decreasing temperature 
in these compounds. Obviously a ferromagnetic exchange in- 
teraction of a localized spin and conduction electrons is the 
first candidate. In fact, later study on dilute Gd and Nd impu- 
rities in some La alloys such as LaA^B, and LaSn^l 3 ] explained 
the resistivity anomaly in terms of the ferromagnetic exchange 
modeltj, and was referred to as the "reverse Kondo effect". 

In dilute Rh alloys, however, the measured susceptibility 
indicates antiferromagnetic exchange interaction in these ma- 
terials in contrast with the behavior of the resistivity. As al- 
ternative interpretation, Fischer found that a strong potential 
scattering caiuchange sign of the Kondo logarithmic term in 
the resistivityQ. The_treatment has been much extended and 
deepened by KondoO, who uses the scattering phase shift S v 
of conduction electrons at the Fermi surface. The strength 
v of the potential scattering is related to the phase shift by 
tan 5 V = —irvpc where p c is the density of states at the Fermi 
surface. The leading logarithmic term in the electric resistiv- 
ity changes sign when \S V \ exceeds 7r/4, defining a critical 
value for the potential scattering as v CT = 1/ (7rp c ). The range 
\v\ > v CI is called reverse Kondo range. In this range the resis- 
tivity decreases with decreasing temperature as a consequence 



of the strong potential scattering. Regarding thermodynamic 
properties, Kondo showed that the effect of ordinary scatter- 
ing is entirely absorbed into an effective exchange interaction 
J = J cos 2 S v . Correspondingly, the low-temperature energy 
scale is characterized by temperature 7k = De 1 ^ 2pcJ ' ) asso- 
ciated with the effective exchange interaction J. 

Since ordinary scattering events are always present in real 
systems, the Kondo problem with strong potential scatter- 
ing might have relevance also in other compounds that show 
Kondo effect. For example, the question of relevance of 
ordinary scattering in URu2Si2 arises, because: (i) recent 
STM experiments have found that the density of electronic 
states shows Fano lineshape, i.e. antiresonance, in the nor- 
mal phase,!] and (ii) in the dilute system U x Thi_ JRu2Si2 the 
resistivity decreases with decreasing temperature.! 

In this paper we study the effect of strong potential scat- 
tering on physical properties of the Kondo impurity. In or- 
der to deal with the Kondo effect beyond the weak cou- 
pling regime, the continuous-Jime quantum Monte Carlo (CT- 
QMC) method is employedMtil. In the CT-QMC simula- 
tion it is most conveniejut-to take the iV-component Coqblin- 
Schrieffer (CS) modetlL 2 ] with potential scattering. The 
Hamiltonian is given by 



n CS [vcs} = X! £kC km C km 
km 

+ [ J fnJm> + V C s5 mm >) <4,C m , (1) 



where and are creation operators of conduction and 
localized electrons, respectively, at the impurity site with 
SU(iV) index m = 1, . . . , N. The constraint J2 m fLfm = 1 
is imposed, which removes the charge degrees of freedom. 
The annihilation operator c m in the Wannier representation is 



related to Ckm by c„ 



-1/2 



J2w c km with N being the 



2 



number of lattice sites. We observe the relation 

^ fmfm' C m' C m = ^ mm ' C m' C m + • 



2S/-s c + -n c , (2) 

JV=2 2 



where X mm / = /^/ m , -S mm >/N are SU(A^) generators, and 
s c and n c are spin and charge density operators of conduction 
electrons at the impurity site. The SU(7V) Kondo Hamiltonian 
Hk [v] with potential scattering v is introduced by the relation 



fics[vcs] = %cb = v C s + J/N], 



(3) 



in view of eq. (0). Some typical cases of the model given by 
Hamiltonian (|J) are 

(i) The conventional CS model with vcs = 0, or v = J/N; 

(ii) The SU(A r ) Kondo model with v = 0, or vcs — 

-J/N; 

(iii) Reverse Kondo range with \v\ > v CT = 1/ (^p c )- 

On the basis of accurate numerical results for strong po- 
tential scattering, we investigate the reverse Kondo range in 
detail. Furthermore, properties are studied by changing the or- 
bital degeneracy TV for the CS model. This paper is organized 
as follows. In Section [n| numerical results for the magnetic 
susceptibility are presented. The characteristics of the impu- 
rity i-matrix are discussed in Section III Furthermore, numer- 



ical results are given for transport properties. Section IV is de- 
voted to discussion of quasi-particle properties including the 
Friedel sum rule and optical theorem. The summary of this 
paper will be given in Section |v[ 



II. MAGNETIC SUSCEPTIBILITY AND UNIVERSALITY 
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FIG. 1: Upper: Temperature dependence of the static susceptibility 
for the SU(iV) Kondo model with potential scattering. Potential scat- 
tering v and exchange J are chosen as v — 0,-0.85 ( J = 0.3), and 
v = -0.48 (J = 0.44) for N = 2, and v = 0.009 (J = 0.075), 
and v = 0.0125 (J = 0.0125) for N = 8. Lower: Kondo tempera- 
ture for several values of potential scattering obtained in simulation. 
The theoretical result Ik = De 1 ^ 2pJ ' > is also shown as dashed line. 



First, let us discuss the behavior of the magnetic suscepti- 
bility for a given value of the orbital degeneracy N. The static 
susceptibility is obtained from the imaginary time data by in- 
tegration as 

X (T) = / dr X (T) = / dT(T T M H (r)M), (4) 
Jo Jo 

where the dipole moment M is given by M = ^ Q TO a/l/ Q 
with coefficients m a chosen as J2 a m an = 0, and the super- 
script H denotes the Heisenberg pictureQ. We use a constant 
density of states for the conduction electrons in the simulation 

as 

p c (e)=p Q(D-\e\), (5) 

where po = 1/ (2D) with D = 1 as a unit of energy. In the 
numerical study, we determine the Kondo temperature from 
the low-temperature static susceptibility as 

T K 1 = X (T 0)/Cjv, (6) 



where Cn is the Curie constant. The critical strength v CI is 
given by v cr = l/(npo) = 0.637. 

Figure |l] shows \{T) of the SU(7V) Kondo (or CS) models 
with several values of the potential scattering v and orbital de- 
generacy N (upper), and 7k for different values of v (lower). 
The result Tk = Da 1 / ^ 2p ^ obtained by Kondoi is also shown 
for comparison. We observe in Fig. [l] that the susceptibility 
shows universal behavior as a function of T/Tk independent 
of the value of the potential scattering. Note that the data with 
v = —0.85 for N = 2 as shown by blue symbols are in the re- 
verse Kondo range with S v being larger than tt/4. Even in this 
case, the temperature evolution of the magnetic susceptibility 
shows the universal behavior. 

Now, let us turn to discuss the properties by changing the 
orbital degeneracy 7Y. For N = 2, the static susceptibil- 
ity decreases monotonically with increasing temperature. On 
the other hand, for large values of the orbital degeneracy like 
N = 8 in Fig. [l], the susceptibility first increases as the tem- 
perature is increased, and then decreases in accordance with 
the free moment behavior \ ~ ^/T for large temperatures. 
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The initial increase can be understood in terms of the den- 
sity of states of quasi-particles. To illustrate the mechanism 
of the increase, let us consider the case of v = in the non- 
interacting Anderson model, which simulates qualitatively the 
quasi-particle density of states. Using the Sommerfeld expan- 



X(T) 



Pf(e) 



df(e) 
de 



de 



Pf(0) 



1 



6 Pf (0) 



0{T A 



(7) 



with the /-electron density of states pf(e), we obtain at low 
temperatures 



A 2 ) 



A 2 



C(T 4 ) , (8) 



where r\ and A are the shift and the width of the resonance 
peak appearing in the density of states Pf(e) at low tempera- 
tures. 

For N = 2 the resonance peak is centered at the Fermi 
energy, which gives rj = 0. Therefore, the coefficient of T 2 
in x{T) is negative, i.e. the susceptibility decreases with in- 
creasing temperatures. Increasing the value of degeneracy N, 
the resonance moves to higher energy above the Fermi level, 
i.e. rj ~ Tk, while its width narrows as A - T K /N.U Thus, 
the coefficient of T 2 in expression ([^) becomes positive, so 
that the susceptibility first shows increasing behavior as the 
temperature is increased. 



III. TRANSPORT COEFFICIENTS 

A. t-matrix and Fano lineshape 

We have already shown in eq. (||) that the Kondo and CS 
Hamiltonians are related to each other through a potential 
scattering term. In the simulation for the CS model, instead 
of the bare Green's function g, another Green's function gcs 
is used that absorbs the potential scattering vqs'- 



gcs = 5/(1-- vcsg), 



(9) 



where g(z) = J2k( z ~ £ k) _1 - Then the simulation gives the 
renormalized Green's function G of conduction electrons. We 
introduce a quantity tcs by the relation 



G — gcs + gcstcsgcs- 



(10) 



On the other hand, the t-matrix t of conduction electrons is 
defined by the relation 



G 



gtg- 



(ii) 



By comparing eqs. (|K)|) and (|TT]), we obtain t from tcs by the 
relation 



t = ucs/(i - v C sg) + £cs/(i - vcsg)' 



(12) 



In the special case of v = vcs + J/2 = 0, we recover eq. (33) 
in Ref. | 

For the moment, we concentrate on the case of N = 2. 
In the CT-QMC simulation, the t-matrix is derived in the 
imaginary-time domain. In order to obtain properties in the 
real energy domain, analytic continuation of the numerical 
data is done by using Pade approximation. Figure || shows 
the energy dependence of the impurity t-matrix for three dif- 
ferent values of the potential scattering at various tempera- 
tures. To simplify the notation, we take the convention in this 
paper that energy including an infinitesimal imarginary part, 
e + id, is simply written as e. The left panel of Fig. ^ with 
the value v = 0.15 corresponds to the CS model, center panel 
with v = to the ordinary single-channel Kondo model, while 
right panel to a strong potential scattering > v cl - = 0.637. 
In the case of the ordinary single-channel Kondo model the 
spectrum is symmetric with respect to the Fermi energy, be- 
cause the model has particle-hole symmetry in this limit. 
Increasing the value of the potential scattering, the Kondo 
peak first moves to higher frequencies above the Fermi en- 
ergy. Finally, for strong potential scattering the spectrum be- 
comes highly asymmetric showing an antiresonance around 
the Fermi level. 

Interpretation of the asymmetric spectrum with large \v\ 
can be provided in terms of the Anderson model, which re- 
produces the CS model (|l]) in the limit of deep local electron 
level £f and large Coulomb repulsion U as compared with 
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hybridization V. Namely we take the limits Sf — > — oo, 
Ef + U —> oo and V 2 po —> oo, keeping the ratio J = 
— 2V 2 po/sf finite. Now we construct the /-electron Green's 
function G f v of the Anderson model in the presence of po- 
tential scattering. Let us first consider the pure case v = 0. 
Introducing the irreducible part F(z), we obtain the Green's 
function 



G f (z) = F(z)[l + V 2 g(z)G f (z)}. 



(13) 



In the presence of potential scattering, the /-electron Green's 
function Gf v (z) satisfies the following relation 



G fv (z) = F v (z)[l + V 2 g v (z)G fv (z)} 



(14) 



where g v (z) = g(z)/[l - vg(z)]. 

The i-matrix for the CS model is given by 



t(z) 



v + V 2 F v (z) 



l-g(z)[v + V 2 F v (z)} 



t v {z) + 



V 2 G fv (z) 
[1 - vg{z)f 



(15) 



where t v = v/(l — vg). It is clear from eq. ( |15[ ) that the 
i-matrix reduces to 

t(z) -> V 2 G f {z) (16) 
in the limit of v = as we expect. We derive from eq. (Iff 



V 2 G fv = {l-vgf{t-t v ), 



(17) 



which is valid for any value of parameters. 

As the simplest case, let us consider the non-interacting 
Anderson model with U — 0. Then we obtain F(e) = 
l/(e + iS — St ) = l/£, which is not affected by potentital 
scatteing. The imaginary part of the /-electron Green's func- 
tion is given by 



i™r (f) — — *po(v + v2 /0 2 

M j ~ V 2 ' 1 + 7r 2 p 2 (v + V 2 /i) 2 ■ 



(18) 



Under the transformation v — > —v, the imaginary part of the 
Green's function given in eq. ( |l8| ) remains the same if we put 
e — Ef^rEf — e. Physically it means that the antiresonance is 
reflected with respect to the energy e f when v changes sign. 
Expression. (|I^) can be put into the standard form of the Fano 
lineshapeo. Introducing the dimensionless parameter q by 
l/q = nvpo, we rearrange the terms as follows: 



V 2 

Im Gj v (e) 



1 



(x + q) 2 



q + l/q x 2 + 1 
where x is the dimensionless energy defined by 



q+ l)t + -q 



(19) 



(20) 



The degree of asymmetry is determined by the parameter q 
that is independent of hybridization. 

Expression ( p"8[ ) describes the characteristics of the simula- 
tion results for t (e) shown in Fig. |[ provided we put e/ ^0. 
We note that V 2 Gf v (z) is not the same as the t-matrix t{z) of 
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FIG. 3: Temperature dependence of the normalized resistivity for the 
Kondo model with potential scattering for N — 2. Potential scatter- 
ing v and exchange J are chosen as v = 0, —0.85 (J = 0.3), and 
v = 0.78, -0.48, -0.78 (J = 0.44). The dashed line corresponds 
to Kondo's result for the resistivity given in Ref . |i| 



the CS model as shown in eq. (|15|). However, the character- 
istic lineshape comes almost from V 2 Gf v (z) since t v (z) and 
1 — vg(z) do not have strong dependence on z. The asym- 
metric spectrum for strong potential scattering has interesting 
consequences with respect to transport properties such as the 
resistivity or thermopower. This problem is discussed in the 
next subsection. 



B. Relaxation time and transport coefficients 

We rely on the Boltzmann equation approachEl to derive 
transport coefficients. Then the relaxation time t(e) is related 
to the i-matrix asLj 



r(e)" 1 = -2Imt(e). 
Let us introduce the integralsQ 



(21) 



(22) 



in terms of which the conductivity er, thermopower S and ther- 
mal conductivity k are expressed as 



<t(T) = Lq, 



S(T) = - 



lit 



(23) 
(24) 

(25) 



The integrals in eq. ( |22[ ) are evaluated numerically with the 
CT-QMC data for the i-matrix. 

Although many theoretical attempts were made to derive 
the transport properties for the Kondo problem analytically in 
the whole temperature range, none of these attempts was suc- 
cessful. However, there are correct results in limiting cases 
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la0 



such as. Hamman's formulaliil for T ^> Tk, Fermi-liquid 
resultslia for T <C TWiOr expressions for large values of the 
orbital degeneracy iVJl3 In the local Fermi-liquid range at low 
temperatures, the resistivity R(T) shows the following tem- 
perature dependence: 



R(T)/R(0) = l-a(T/T K y 



(26) 



where we have used the relation R(T) /R(Q) = cr(0) /cr(T), 
and a is a numerical coefficient. For large values of the orbital 
degeneracy_/V, the l/N expansion gives the coefficient a in 
eq. © asO 



a = 7r 



1 



3N 



(27) 



This limiting result gives checkpoint of our numerical calcu- 
lations. 



C. Behavior under varying potential scattering 

Figure |] shows the temperature dependence of the normal- 
ized electric resistivity across the Kondo temperature. For 
\v\ < v CI = 0.637, the resistivity follows universal behav- 
ior as a function of T/Tk. As the potential scattering in- 
creases beyond the critical value, the resistivity still follows 
the universal behavior in the Fermi-liquid range T <C 7k- 
As temperature increases, the resistivity starts to deviate from 
the universal curve around the Kondo temperature T ~ Tk, 
and shows increasing behavior as the temperature is further 
increased. In the temperature range T ^> 7k, we find that 
the resistivity for strong potential scattering can be described 
well with Kondo's formulas shown by dashed line in Fig. 
However, Kondo's formula cannot describe the properties for 
T < Tk- Based on the numerical results for T <C 7k, the 
coefficient a in eq. (^) appears to be independent of v. 

The upper panel of Fig. ^ shows the temperature depen- 
dence of normalized thermal conductivity for N = 2 under 
varying the potential scattering. The thermal conductivity also 
follows the universal behavior even for large values of the po- 
tential scattering in the temperature range T < Tk. Namely 
we obtain 



k(T) . (k{T) 



T 



-I 



T 



1 + 7 



T 
Tk 



(28) 



with 7 being a numerical constant independent of v. Increas- 
ing further the temperature, thermal conductivity with large 
potential scattering highly deviates from the universal behav- 
ior. Namely, it decreases with increasing temperature for 

T>T K . 

The lower panel of Fig. [| shows the temperature depen- 
dence of thermopower for N = 2 under varying the potential 
scattering. The asymmetry of the impurity i-matrix is most 
reflected in the behavior of thermopower. This is clear if we 
regard explicitly the expression for the integral L\ given in 
eq. (§2h: 



L 



f'(e)e 
lmt(e)' 



(29) 
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FIG. 4: Temperature dependence of normalized thermal conductiv- 
ity (upper) and thermopower (Jower) for the Kondo model with po- 
tential scattering for N — 2. Potential scattering v and exchange J 
are chosen as v = ( J = 0.3), and v = 0.78, -0.48, -0.78 (J = 
0.44). The dashed line in the right panel corresponds to the Fermi- 
liquid behavior S(T) ~ T. 



Because of the factor e in the numerator, the thermopower 
measures the asymmetry in the energy dependence of Im t(e). 
Since the spectra is completely symmetric in the case of the 
ordinary Kondo model (v = 0), the thermopower vanishes in 
this case as we obtain in the simulation (see Fig. ||). On the 
other hand, the thermopower acquires strong temperature de- 
pendence when the potential scattering term is increased from 
v = 0. The thermopower in the simulation shows Fermi- 
liquid property 



(30) 



for T <C Tk- The coefficient /3 is negative for v < 0, while 
it is positive for v > 0. The different sign of (3 can be simply 
understood if we recall that the asymmetry in the lineshape 
of Imt(e) around the Fermi level is reversed against the sign 
change v — > —v. Thus, the integrals given in eq. ( |22| ) have 
the properties Lq(L2) — > Lq(L2) and L\ — > —L\ under the 
sign change. Namely, the electric resistivity R(T) and thermal 
conductivity k(T) remain the same under v —v, while 
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FIG. 5: Upper: Temperature dependence of the normalized resis- 
tivity for the CS model (where v = J/N) for orbital degeneracies 
TV = 8 with J = 0.075, 0.1 and TV = 50 with J = 0.0115. 
Lower: Orbital degeneracy TV-dependence of the coefficient a of the 
T 2 term in the low- temperature resistivity defined in eq. (^). The 
inset shows the scaling behavior of the resistivity for different orbital 
degeneracies including a. 



the thermopower changes sign as S(T) — > —S(T). We have 
indeed obtained these behaviors in the simulation as it can be 
seen in Figs. [3] and 



D. Behavior under varying orbital degeneracy TV 

Let us now study the effect of orbital degeneracy. The upper 
panel of Fig. |] shows the temperature dependence of the nor- 
malized electric resistivity for large values of the orbital de- 
generacy for the CS model with v = J/N across the Kondo 
temperature. We observe again the universal behavior for a 
given value of the orbital degeneracy TV. In contrast to the 
behavior of the magnetic susceptibility shown in Fig. |IJ the 
resistivity decreases monotonically as the temperature is in- 
creased even for large TV. In order to understand this feature, 
we assume that the i-matrix at low temperature is determined 
by the quasi-particle density of states, which is approximately 



given by the effective Anderson model. Namely we assume 

t{z) = V 2 G f {z) 1 (31) 

where V is the effective hybridization and G / is the Green's 
function of the local electron in the effective Anderson model. 
Then the Sommerfeld expansion of the conductivity leads to 



v(T) 




where 



P/(e) 



Mo) 

(32) 



(33) 



We obtain the resistivity R(T) = cr(T)" 1 from eq. (||). Us- 
ing the quasi-particle density of states for the non-interacting 
Anderson model with v = 0, we obtain 



R{T) ~ pf(ft) 



A 2 



0{T 4 



,(34) 



where r\ and A are the shift and the width of the resonance 
peak appearing in Pf(s) at low temperatures. Irrespective of 
the magnitude of parameters r\ and A, the coefficient of the T 2 
term in the low-T resistivity is always negative. Hence R(T) 
given by eq.(Q) decreases as temperature increases for any 
value of TV. Namely, the quasi-particle picture is consistent 
with the monotonous change obtained in the simulation. 

The resistivity obtained numerically has T 2 temperature de- 
pendence at low temperatures, which is expressed generally in 
eq. (|^). The lower panel of Fig. |] shows the coefficient a of 
the T^term for several values of the orbital degeneracwJV ob- 
tained in the simulation. The result of 1/TV expansion^ given 
in eq. (|7|) is also shown in the figure. We find that the nu- 
merical data coincide with the 1/TV expansion result for large 
values of the degeneracy TV. Thus, our CT-QMC simulation 
has produced accurate numerical results for large values of 
orbital degeneracies TV — > oo, which might be difficult in the 
case of other numerical techniques. For small values of de- 
generacy TV, the coefficient a shows linear- TV dependence in 
our simulation. 

Now we fit the simulated TV-dependence of a by a rational 
function as a(N) = F x (0 (TV)/F 2 M (TV), where = c k0 + 
CfeiiV+ ...CkiN 1 . We find that the minimal function which can 
give a good fit to the numerical data has the form 



a(N) = 



Cio + ciiTV + ci 2 TV 2 
1 + c 2 iN + C22N 2 



(35) 



The limiting cases TV — > and TV — > 00 are obtained from the 
formula (Ba) as 



a(N -> 0) ~ C10 + (en - ci c 2 i)TV 



(36) 
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FIG. 6: Temperature dependence of normalized thermal conduc- 
tivity (upper) and thermopower (lower) for the CS model (where 
v = J/N) for orbital degeneracy N = 8. Potential scattering terms 
are chosen as v = 0.009 (J = 0.075) and v = 0.0125 J = 0.1). 
The dashed line corresponds to the Fermi-liquid behavior given in 
eqs. (5(3) and (§§). 
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(37) 



Choosing the coefficients in eq. (g5J) as c 10 = 0.3; c xl = 0.53; 
c 12 = 0.17; C21 = 0.1; c 2 2 = 0.0 1 73, the numerical data can, 
be fitted well (see Fig. ||). The result of l/N expansionli-1 
given in eq. ( |27| ) is also reproduced in the large- range. In 
the inset of right part of Fig. || we plot the normalized resistiv- 
ity for different values of orbital degeneracy TV as a function 
of a(A r ) 1 / 2 T/TK- The results show the scaling behavior of 
the resistivity in the Fermi-liquid range. For T > 7k, the 
scaling property breaks down. 

In Fig. [| thermal conductivity k and thermopower S are 
shown for the CS model with orbital degeneracy N = 8. 
The low-T behaviors of k(T) and S(T) are consistent with 
the Fermi-liquid result given in eqs. ( |2q ) and (p0|). As the 
temperature is further increased, the thermopower S(T) has a 
peak, while the thermal conductivity k(T) monotonously in- 
creases. Both quantities show universal behavior as a function 
of T/T K . 



Note that the coefficient j3 for S(T) has a positive sign for 
large N. This behavior is explained as follows. In a manner 
similar to eq. (f32|), L\ is given by 



-00 P/(e) 



3A 



(38) 



Using the result for Lq given by eq. (p4[), we obtain S(T) in 
the lowest order of T as 



S{T) 



2tt 2 



V 



>r 



r T = (3T. 



(39) 



Since r\ ~ Tk > for large TV, we obtain f3 > 0. On the other 
hand, we obtain i] = for the symmetric Anderson model 
with TV = 2. In this case, the sign of j3 depends on the sign of 
v as it can be observed in Fig. H 



IV. QUASI-PARTICLE PROPERTIES 
A. Friedel sum rule 

At temperatures T -C Tk, the conduction electrons screen 
the magnetic impurity and they together form a local singlet. 
In this range the ground state is a local Fermi-liquid. The 
Friedel sum rule (FSR) relates the phase shift for scattering of 
the conduction electrons by the impurity to its charge. In the 
case of the CS model, the /-electron Green's function cannot 
be defined because there is no charge degrees of freedom in 
this localized model since it is eliminated. Instead, the im- 
purity t-matrix is used to describe the effect of exchange and 
potential scatterings. In Section Q we have related the t- 
matrix to the Green's function Gf v {z) of localized electrons 
with potential scattering v. It is the quantity Gf v (z) that is 
expected to keep the FRS in the presence of v. 

The FSR readsB 



V 2 G f (0) = -^sin 2 (£) , 
npo \ TV / 



(40) 



since the occupation number is unity in the CS and Kondo 
models. 

Figure ^ shows the Green's function Gf v (0) obtained by 
simulation at finite temperatures. For small potential scatter- 
ing the simulation results show good agreement with the ex- 
pectation given in eq. (ftfj|). As the value of v is increased, 
Re V 2 Gf v (0) is still close to zero, but -Im V 2 G fv (Q) highly 
deviates from the theoretical result. The reason of this devia- 
tion is the following. The theoretical result given in eq. ( |40[ ) 
is realized at T = 0, which is almost realized for T Tk 
in the simulation. As the value of the potential scattering 
is increased, however, the corresponding Kondo temperature 

Tk = Dc 1 /( 2pJ * > rapidfvL decreases since the effective cou- 
pling J decreases rapidlyQ Therefore, we have to go at lower 
and lower temperatures in the simulation to achieve the condi- 
tion T ^Tk, which is not fulfilled in Fig. [7]for large values of 
v since the results are obtained for a fixed temperature value 
/3 = 1/T = 1000. In principle it is possible to go at lower 
temperatures in the simulation, but it becomes computation- 
ally harder. 
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FIG. 7: Imaginary and real parts of the Green's function Gf v ex- 
pressed in eq. ( Jiv| ) at the Fermi level as a function of potential scatter- 
ing v together with the analytical result for — Im V G f v (0) obtained 
from the FSR (dashed line). The figure also shows the inelastic scat- 
tering cross section (Tmci as blue triangles. The numerical data are 
obtained at temperature /3 = 1/T = 1000. 



B. Optical theorem 

The optical theorem is less restrictive than the FSR since 
the former does not require the Fermi liquid ground state. Op- 
tical theorem is related to the unitarity of the 5* matrix,E3 and it 
follows when the scattering of the conduction electrons from 
the impurity is totally elastic at the Fermi level. Optical the- 
orem was originally formulated for problems of scattering of 
a single-particle. When the scattering event happens without 
energy loss, i.e. it is totally elastic, there is a relation between 
the square and the imaginary part of the t-raatrix. To express 
this relation, the S-matrix is decomposed ascll 



S = 1 + iT, 



(41) 



where we write the matrix element of the t-matrix T as 
{n\T\n') = 2TrS(e n — e n ')(n\t\n'). Here, a state \n) repre- 
sents a single-particle state with momentum k and spin tr as 
\n) = | kg). After some manipulations we obtain from rela- 
tion (0) that 



(k\SS*\k) - 1 = 2tt 



2nJ2S{e k ~e n )\(n\t\k)\' 



2lm(k\t\k}] 



(42) 



The first term of eq. (^2j) in the right-hand side is related to 
the elastic scattering cross section cr e i, while-lhe second term 
to the total scattering cross section cr t otai asElo 

CT el = -27r^«5( £fc -e„)|(n|t|fc)| 2 , (43) 

71 

atotai = 2lm(k\t\k). (44) 

The inelastic scattering cross section <7i no i is the difference of 

CTtotai and <j cl as 



For scattering only in the s channel and assuming 
servation, the inelastic cross section is expressed as 



=iM = OH 

= 27TPQ 



1)/(2tt) = 2Imi(w) + 27r ( o |t(u;)| 2 
— Imt(w) + \t{oj)\ 2 , (46) 

7T/00 



where s and t are eigenvalues of the S'-matrix and i-matrix, 
respectively. 

The eigenvalues s of the S'-matrix lie within the complex 
unit circle. The scattering is completely elastic, i.e. <7i nc i = 
0, when the unitary condition S S' = 1 is satisfied, which 



means s 



1. In this case we have the relation 



-Imt 



(47) 



from eq. (Q. 

The problem is formulated so far for single-particle scat- 
tering, but it is more general and can be applied for many- 
particle problems as well such as the single-channel Kondo 
model. It is the most easily understood in the limit of e — > oo, 
when the magnetic impurity is completely decoupled from the 
conduction electrons. In this case the conduction electrons 
scatter without energy loss, and the optical theorem is held. 
Although the t-matrix is complicated and contains many scat- 
tering events at low temperatures, strictly in the limit of e — > 
the optical theorem holds again. If we express the t-matrix as 

t = \t\e* e , (48) 

where 6 is the phase of the ^-matrix t. Equation (^7|) gives 



Imt(e) 



sin 



7T/O0 



(49) 



for the Kondo model with energy not only e = but also 
e — > oo. Namely, relation ([l9]) is satisfied in the case when 
the scattering is totally elastic. In the case of ordinary single- 
channel Kondo model (v = 0), 9 = —ir/2 at the Fermi energy, 
so from eq. (B9h we recover the FSR 



lmt(0) 



1 



(50) 



Cinel = CTtotai CT e i- 



(45) 



Let us discuss the optical theorem in the context of our nu- 
merical data. Figure [7] shows the obtained inelastic scattering 
cross section a- lne \ as a function of potential scattering. Note 
that CTinei shows non-monotonous behavior as a function of v. 
Namely, aj ne i is almost zero for small values of the potential 
scattering, but becomes non-zero as \v\ is increased. With fur- 
ther increase of \v\, however, er inc i approaches to zero again. 
Since Tk decreases as the potential scattering increases, in the 
large-w range the condition T 3> Tk is satisfied at T = 10~ 3 
used in the simulation. On the other hand, in the small-w 
range we have the condition T <C Tk with the same value: 
T = 10~ 3 . It is confirmed in the simulation that the opti- 
cal theorem holds both at T < T K and T > T K , but not 
for T ~ Tk- This happens because the ranges T <C Tk and 
T > T K correspond to the limits£/T K < 1 and e/T K > 1, 
respectively, where expression ( f49| ) is satisfied as explained 
above. 
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V. SUMMARY 

In this paper we have studied Kondo impurity models with 
potential scattering and orbital degeneracy by using CT-QMC 
numerical technique. We have derived accurate numerical re- 
sults for the impurity i-matrix, thermal, and transport proper- 
ties in a wide temperature range across the Kondo tempera- 
ture Tk- Properties in the reverse Kondo range has been in- 
vestigated in detail. The results shown in this paper are nu- 
merically exact since CT-QMC does not use any approxima- 
tion. We have explicitly demonstrated that CT-QMC simula- 
tion technique gives numerically exact results for large val- 
ues of the orbital degeneracy N, which might be difficult to 
achieve in the case of other numerical techniques. 

For large values of the potential scattering, non-trivial 
physics appears even in the impurity problem. Namely, the 
resistivity shows anomalous increase with increasing temper- 
ature in contrast to the ordinary Kondo effect. This unusual 
behavior is caused by an antiresonance developing around the 
Fermi energy in the quasi-particle density of states as the value 
of the potential scattering is increased. This antiresonance 
does not influence the universal behavior of the magnetic sus- 
ceptibility. However, the sign of the Kondo logarithmic term 
changes in the resistivity when the potential scattering exceeds 
a critical value, i.e. in the reverse Kondo range, which causes 
the resistivity decrease with decreasing temperature. 

We have studied the effect of strong potential scattering on 
thermal and transport properties of the Kondo impurity, and 
obtained that 

(i) the magnetic susceptibility follows the universal tempera- 
ture dependence even with strong potential scattering; 

(ii) the resistivity also follows the universal temperature de- 
pendence for small values of the potential scattering; 

(iii) when the potential scattering exceeds a critical value, the 
resistivity still shows universal behavior in the Fermi-liquid 
range, but starts to deviate from the universal curve around 
the Kondo temperature, and increases as the temperature is 
further increased; (iv) the obtained temperature dependence of 
the resistivity for T ^> 7k in the reverse Kondo range agrees 
quantitatively with Kondo's theoretical result; 



(v) the thermal conductivity also shows universal behavior in 
the Fermi-liquid range, but highly deviates from the universal 
curve for T 3> 7k in the reverse Kondo range; 

(vi) the asymmetry of the t-matrix developing with increasing 
value of the potential scattering is most reflected in the tem- 
perature dependence of the thermopower; 

(vii) the sign of the thermopower depends on the sign of the 
potential scattering. 

In addition to the study of thermal and transport properties, 
we have discussed the Friedel sum rule and optical theorem as 
well. We have shown that the t-matrix of the Kondo model in 
the presence of potential scattering is not the relevant quantity 
for the Friedel sum rule. Instead, the Friedel sum rule is sat- 
isfied with a proper limit of the /-electron Green's function. 
We have demonstrated that optical theorem is less restrictive 
than the Friedel sum rule, because the former holds not only 
in the Fermi-liquid range, but for large energies as well. 

Finally we mention an interesting question whether the be- 
havior found for strong potential scattering has relevance in 
real systems. Note that recent STM experiments on URu2Si2 
have found that the density of states shows Fano lineshape in 
the normal phase, and in the dilute system U T Thi_ x Ru2Si2 
the resistivity decreases with decreasing temperature. Since 
some important aspect of the U ion with non-Kramers con- 
figuration 5/ 2 may not be described by a localized spin with 
S = 1/2, account of the strong potential scattering in more 
realistic models is desirable. We hope that our results in this 
paper will stimulate further study concerning the Fano line- 
shape and other aspects, which reflect interplay of the Kondo 
effect and potential scattering. 
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